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1. Introduction 

In a series of recent papers p El El m Karkowski, Mach and Malec investigated 
spherically symmetric Bondi-type accretion of perfect fluids in Schwarzschild-(anti-)de 
Sitter spacetimes. In P and P a more general case of self gravitating fluid is also 
studied. The primary motivation of these works was cosmological; the authors were 
mainly trying to determine the way in which the presence of the cosmological constant 
influences the possible accretion rate. 

In this paper I go a step further and consider relativistic Bondi-type accretion 
[22| E3] in the Reissner-Nordstrom-(anti-)de Sitter spacetime. The motivation behind 
this choice comes from the fact that the structure of Reissner-Nordstrom metrics 
resembles certain features characteristic for the Kerr solution that could influence 
accretion p. They share similar horizon structure p. Also, the Penrose diagram for 
a generic case of Reissner-Nordstrom-de Sitter spacetime and the Penrose diagram for 
Kerr-de Sitter universe are similar P . These similarities in causal and horizon structures 
let one suppose that accretion solutions in both spacetimes may be similar. In fact it is 
common to treat the Reissner-Nordstrom solutions as a toy model for astrophysical black 
holes PPPP. Finding general, radially dominated, stationary accretion flows on the 
Kerr metric seems to be intractable at present (cf. [T0llTl| [T2]b with the exception of the 
ultra-relativistic potential flows [lEllIl]. Investigating spherically symmetric solutions 
in the Reissner-Nordstrom-(anti-)de Sitter spacetime seems to be a way to avoid (at 
least temporarily) the main mathematical difficulties of dealing with metrics that are 
not spherically symmetric. I believe, however, that steady Bondi-type accretion in the 
Kerr spacetime and in Reissner-Nordstrom-type spacetimes could share some features 
of common behaviour, at least qualitatively. 

Reissner-Nordstrom-de Sitter black holes can also be considered in the context of 
the anti-de Sitter/conformal held theory correspondence [T^ 116] . Thermodynamical 
properties of the Reissner-Nordstrom-(anti-)de Sitter black holes are a subject of 
numerous investigations [HI [ig EH Eo]. 

Very recently Chaverra and Sarbach presented another analysis of Bondi-type 
accretion hows on essentially arbitrary, spherically symmetric metrics [2l] . The solutions 
discussed in this paper fall in their class. Of course, by specializing to a narrower set 
of Reissner-Nordstrom-(anti-)de Sitter spacetimes I was able to obtain much more 
detailed information about corresponding accretion solutions. 

The order of this paper is as follows. In Section 2 I introduce the Reissner- 
Nordstrom-(anti-)de Sitter metric and specify the coordinate system. Section 3 
introduces general equations of the how. I also dehne the notion of the sonic point there. 
Section 4 uses results of the previous section in case of diherent types of isothermal test 
huids. I show, that for some of them it is possible to obtain analytic solutions. I 
plot some of them in a phase diagram of the velocity vs. radius. It shows, that for 
subrelativistic huid and negative cosmological constant it is possible to obtain closed, 
binoculars-like trajectories. Section 5 yields results for polytropic test huids. Also for 
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them there exist closed trajectories. Section 6 summerizes the paper and compares the 
obtained results. In the Appendix I prove the relation between locations of the horizons 
and the sonic points in case of the radiative fluid in Reissner-Nordstrom-(anti-)de Sitter 
spacetime. 


2. Reissner—Nordstrom—(anti—)de Sitter spacetime 


In polar coordinates {t,r,6,(j)) the metric of the Reissner-Nordstrom-(anti-)de Sitter 
spacetime can be written as [T9] : 


= 


2 m 

o 


A 


r r^ 3 


dt^ + 


( 1 - 


2rn \ 0^ _ A^2 

r ' r'^ 3 ' 


dr^ + 


+r\de^ + sm\e)d(j)^), ( 1 ) 

where m and Q correspond to the mass and the charge of a black hole, and A is a 
cosmological constant. This coordinate system is singular at the solutions of the equation 


2m A 2 „ 

-H ^ = 0. 


( 2 ) 


This equation has at most three real, positive roots, corresponding to the Cauchy 
horizon, the event horizon and the cosmological horizon. The last one exists only for de 
Sitter spacetimes (A > 0). In order to get generic case where these three horizons exist 
(or two in case of A < 0) one imposes constraints on the values of m, Q and A. For 
A > 0 there is either A = 0, < 1/(4A) or A > 0, < 1/(4A), and for A < 0 there 

is A < 0, where 


A = (sim^A + Q2(3 + 4g2^) _ 9^2^'^ ^ 12 Q‘^A)') . 


( 3 ) 


I can get rid of the coordinate system singularities at the black hole horizons by 
introducing coordinates with time t' dehned as 


dt' = dt 


2m _ of. I A 2 
r _ ' 3 

\ _ 2m Of. _ A^2 
r 3 


dr. 


( 4 ) 


This yields 

2 f 2 m A 2\ , /2 f 2 ’ 22 - ^ 2\ 7 / 7 

ds^ = -{l -+ ^ - dt'^ + 2 -^ dt'dr + 

\ r r^ 3 J \r r^ 3 ) 

+ — ^ + — dr^ + r‘^{d6‘^ + sin^(6*)d0^). 

\ r r^ 3 J 

The determinant of the metric ([5]) is 5 ^ = —sin^ 9. 


( 5 ) 


3. Flow in Reissner—Nordstrom—(anti—)de Sitter spacetime 

In this article I consider a perfect fluid characterized by the energy-stress tensor 
= (e + p)u^^u’^ + 


( 6 ) 
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where e denotes the energy density, p is the pressure and is the four-velocity of 
the fluid. The Bondi-type (or Michel-type |22l [23]) accretion is a steady, spherically 
symmetrical flow. In both coordinate systems, all quantities should be functions of 
radius only, and = 0. From the normalization of the four-velocity vector I get 


u^' = 


f 2Tn _ : A„2 

r + 3 




_ 2 m I 




A 2 

3 ' 


Uf = -\ I 


2m ^ Q‘^ 


A 


-r- 




(7) 


( 8 ) 


r 3 

One can check, that in both coordinate sytems Ut and Uf has exactly the same form. I 
will describe the motion of fluid using two conservation equations: 


VApaO = 0, (9) 

V„((e+p)ti'‘!i‘' + p9'“') = 0. (10) 

where p denotes the baryonic (rest-mass) density. Equation (I9|) can be rewritten as 

V^(pm'') = -^a^(\/^pn^) = ^dr{r‘^pu'’) = 0. (11) 

Let me introduce a specihc enthalpy h = (e -|-p)/p. Assuming that the flow is 
smooth, equation fITOl) can be expressed as 

= V^ihpu^u^) + g^'^d^p = pu^V^ihv!^) + hu^V ^{pu^) + g^'^d^p = 

= pu^^V,{hu'^)+g^'^d,p = Q, ( 12 ) 

where I used equation (I9|). In the following I will consider isentropic flows, for which 
dh = dp/p, and in consequence 

u^V ^{hu^) + dyh = u^d^{hui/) — T/,^huxu^ + d^h = 0 . ( 13 ) 


Taking v = t' one may hnd out, that the terms containing Christoffel symbols cancel. 
This yields 


drihuf) = 0. (14) 

Note that equations (I9|), fflOl) may be rewritten as equations ffTTll . fflTl) . Let me point 
out, that in both of the considered coordinate systems these equations have the same 
form, and consequently the obtained functions Ut{r) {ut'{r)) and ^^{r) are the same in 
both systems. For this reason I will not differentiate between Ut and Ut' in what follows. 
By integrating equations ffTIl) . ffTTll . and substituting Uf with equation (IH|) one gets 

= const, (15) 

hJl- — + %--r^ + {u^f = const. (16) 

V r r^ 3 

These equations constitute the starting point of the remaining part of this work. 
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Let a denote the local speed of sound. A location in which the four-velocity of the fluid 
satisfies = (pf /utY will be called a sonic point (because of the spherical symmetry of 
the flow, it is actually a sphere). I will consider barotropic equations of state {h = h{p)), 
so 


dh 

— = a —. 
ti p 


(17) 


Taking into account equation (1T7|) . differentiating equations flT^ and (fT6|) yields 


r\ 2 


dr In u"' = 

r[utY 


9/ ^ ^ 9^ 


For a sonic point (if only dr Inn'’ < oo) one has 


m 


2aYut*) h 


Q 2 A 


+ = 0. 


(18) 


(19) 


From now on quantities with an asterisk will refer to the sonic point. Equation ([19]) can 
be rewritten as 


r\2 


(<) 


2n 2rl 6 * *\ 2r, 2r2 2 *) ' 


( 20 ) 


Equations ffTHl) . flTOll and ffTTll have to be supplemented with suitable boundary 
conditions. In the remaining part of the paper the values of the density and the speed 
of sound at the boundary of the cloud are denoted by poo and Ooc- The radius of the 
cloud is denoted by r^o. It can be either finite or infinite depending on the context. 


3.2. Hamiltonian dynamical system 


Differentiating equations ffTHl) and flTOl) with respect to r, and using relation equation 
flT7|l leads to 


du^ 

dr 


2u^ 1 

\ r 

- 

1 - 1 

(m 

\ 2r 2r2 

-y) 

r 

(m^)^ - a2 1 

_ 2m ^ . 

1 r 

- |r2 (■u’’)^1 

1 


One may introduce a parameter /, such that 


dr 

{u^ 

^2 2 
) — a 

fl- 

2m 
-h 


2 1 f r'N 2 

—r + (u 

dl 



V 

r 

J.2 

3 ^ ^ 

dnS 



2m 


A 2 


= 2u' 

1 — 

— 

+ ^ 

- 

+ K) ) - 

dl 


[ V 

r 

^2 

3 



k! _ Ird 

2r2 6 j 


( 21 ) 


( 22 ) 

( 23 ) 


These equations constitute an autonomous, hamiltonian two-dimensional dynamical 
system. Its orbits are composed of the solutions of equations flT^ and flT6l) . The 
parameter I is arbitrary, so the reparametrization of this system shall not change its 
orbits. From equation fl20l) it is clear, that sonic points are the critical points of this 
dynamical system. 










Bondi-type accretion in the Reissner-Nordstrdm-(anti-)de Sitter spacetime 


6 


One may treat left-hand side of the equation fll6p as a hamiltonian H of this system. 
Equation may be recreated from this hamiltonian as 
dr dH 


dl 

du^ 


du''' 

dH 


(24) 

(25) 


dl dr ’ 

although parameters I and I differ. 

Plots in this paper will be given in {H'/utY vs. r variables, as they present 
information in more readable way. One may consider the above dynamical system 
in H jut and r variables. Then it will be described by the hamiltonian 


H = h 


A 

Then equations analogous to equations 


1 _ ^ ^ _ A^2 

r _ rz _ o 


(26) 


and 


yield 




A 2 
3 


dr 


te)' 




A „2 


(27) 


This dynamical system has an additional critical point, which is located at the horizon. 

4. Accretion of isothermal test fluids 

In this chapter I consider the equation of state in the form p = ke where k is constant. 
In this case = dp/de, so = k. Moreover, equation flTTl) yields 


de e + p{e) 

dp p 

By integrating this term from r^o to any point inside ball of fluid I get 


de' 




g \ 1 /(fc+i) 


(28) 


(29) 


where the last equality assumes the equation of state p = ke. Combining equations fl28ll 
and fl29|) yields the following expression for the enthalpy 

(30) 


_ {k + l)eoo 

pk+i P ■ 
roo 

Equation flTOll can be now written as 


k 2m Q2 A ^ , .2 . 

p \ 1 -1- 5 - r^ [u^] = const. 

V r r^ 3 


( 31 ) 


Together with equation flTHll this leads to 
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In the following I will be especially interested in the flows that pass through 
the sonic point, although I also describe non-transonic solutions. Transonic flows 
are interesting for many reasons. They are often considered in a context of spherical 
accretion m [211221 [23]. Their physical meaning comes from the fact that they maximize 
accretion rate [221EU and are global in some cases | 22 |. This second property also holds 
for some types of spacetimes and matter considered in this paper. 

For = k equation f[20l) can be written as 


m 

2r, 


91 

2rl 


A 


rt = k [I — -— 


3m A 


j- 7 -^ 

2r2 2 * 


Equation 


6 * V 2r, 

can be solved numerically for any k. 


( 33 ) 


On the other hand for k = 


1 


3 2 111 
’ 4 ’ 3 ’ 2 ’ 3 ’ 4 


it reduces to at most four degree polynomial equations of u^{r) and can 
be solved analitically. In this paper I will consider four cases: k = 1 (ultra-stiff fluid), 
fc = I (ultrarelativistic fluid), k = ^ (radiation fluid), and k = \ (subrelativistic fluid). 

4 . 1 . Solution for k = 1 

A fluid obeying the equation of state p = ke for k = 1 is called the ultra-stiff fluid. In 
this case equation (j3^ implies that 

1 _ ^ni j_ _ A7.2 

' • (34) 

There exist two solutions u'^{r) differing only in sign. The expression for r* is much 
more complicated. Equation ([33D is identical to the expression for the locations of the 
horizons ([2|) so the sonic point has to be located on the horizon. The solution can be 
obtained analytically, though the resulting expressions are lengthy. Because of that, I 
will not provide explicit solutions. 

Knowing r* allows one to compute ul from equation ([20|). These two numbers allow 
one to obtain the constant C in equation (IMj) : that results in the explicit form of the 
function u^{r). Finally it is possible to get Ut{r) from equation (0). Sample plots of 
{u‘^/ut)‘^{r) for this and other cases are shown in £gures[3l [H 


4-2. Solution for k = ^ 

In this case, describing ultrarelativistic fluid, equation f[55]l results in 
A . n 5m 


4 2 I 

5 ^.-^. + —r.- —= 0 , 


( 35 ) 


It is also quartic, similarly to the k = 1 case, therefore I will not provide the explicit 
solution. On the other hand equation f[3^ becomes a quadratic equation with the 
solution of the form 


= -lor^ ± 1. 

2 2 \ 


2m 

r 


(72^4 _ 4 M _ ZLLl 4- 


r.2 


A^ 

s’ 


(36) 


Again one can obtain r* and from equations f[35|) and f[20|) . and then use them to 
calculate C from equation ([5BD . repeating procedure from the previous paragraph. 
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With two solutions (15^ . for a generic set of parameters I get two possible functions 
{u^/ut)'^ describing fluids passing through the sonic point. This corresponds to standard 
nonrelativistic accretion considered by Bondi in [22]. The standard interpretation is that 
one of the mentioned functions describes gas accreting onto black hole, while the second 
one is connected with so-called stellar wind. Both solutions are shown in hgures E] 01 
Small changes of the constant C lead to solutions characterized by different values 
of entropy. An example of solutions generated for other values of C, and therefore not 
passing through the sonic point, is presented in figure 01 

(l/lUtf 



Figure 1. Solutions obtained for isothermal equations of state p = e/2 with m = 1, 
Q = 95/100 and A = 1/25 for C constant (from equation (l36l) l equal to 0.146, 0.246 
(transonic solution), and 0.346. Dashed vertical lines denote locations of the horizons. 


4-3. Solution for k = ^ 

The case /c = | is especially interesting, as it describes radiation fluid. Moreover, for 
/c = I equation fl3^ can be reduced to the quadratic equation 

rl - 3mr, + 2Q‘^ = 0. (37) 

The critical points are located at 

3m ± \/9m^ - 8 Q‘^ 

r* = -^-• (38) 

Please note, that the location of the critical point is independent of the value of the 
cosmological constant. The same effect can be observed for the Schwarzschild-(anti-)de 
Sitter spacetime [2] . As it can be seen in hgures El 0] one of these two points is a saddle 
point, when the second one seems to be a center point. 
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The relation between radial component of the fonr-velocity and r can be 
described as 

2 m A 2 / r\2\ ^ 4 / r\2 (‘iVO 

1-^ ^ + (« ) =Cr {u) . (39) 

r r^ 6 J 

As before one can calculate 

^ ,^2 _ ±3mV9m^ - 8 Q'^ (1 - AAQ^) - 9nR ( 4 AQ 2 + i) + 4 Q 2 ( 4 ^g 2 ^ 3 ) 

- 4gg2 • ^ 

Inserting equations (155]) and (HU]) into equation (1551) allows one to hnd the value of the 
constant C depending on m, A, and Q. Varying the value of C leads to results similar 
to the ones for the fc = | case. 

There exists a correlation between the locations of the critical points for k = ^ and 
black hole’s horizons, which is shown in hgure |2l It can be proven, that if there exist 
three horizons, one of them must be above both critical points, one of them must be 
between the critical points, and one must be below both critical points. For A < 0 , when 
there exist two horizons, one is below both critical points and one is located between 
them. The proof can be found in the Appendix. This property implies, that the horizons 
of two types can overlap each other only in the critical points. 


4■4- Solution for k = \ 


When k = equations 
A 


and 


7m 


5Q' 


yield 

= 0 , 


—-—r* + 

6 2 


' 2 m A 2 

r r^ 6 


= Cr^u 


2„,r 


(41) 

(42) 


By repeating the procedure from previous paragraphs one can obtain the function 
(p//utY{r). Again varying the constant C in equation (1521) results with solutions that 
are similar to those obtained in the previous case. 

The solutions are showed in hgures 15], H] One can observe, that in this case it is 
possible to get a closed trajectory on a plot {u'^/utY vs. r, as it is depicted in hgure 01 
Further discussion can be found at the end of the next section. 


5. Accretion of polytropic test fluids 


In this section I assume a polytropic equation of state of the form p = Kp^, where K 
and F are constants. For this equation of state 
F - 1 


h = 


V-I-a? 
Together with equation (1T6|) it gives 


F - 1 - a2 


= const. 


(43) 


(44) 
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r 



Figure 2. Locations of the horizons and the critical points depending on the 
cosmological constant A and the black hole’s charge Q for the equation of state in 
a form p = e/3. Solid, red lines denote Cauchy horizons, dashed, blue lines denote 
event horizons, dotted, black lines denote cosmological horizons and a thick, black line 
denotes critical points. Different lines correspond to cases with different values of the 
cosmological constant A (changing from -16/90 to 1/5). Location of the critical points 
does not depend on the cosmological constant. 





Figure 3. Transonic solutions obtained for the isothermal equation of state p = ke 
with m = 1, Q = 95/100 and A = 1/100. Dashed vertical lines denote locations of the 
horizons. 
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Figure 4. Transonic solutions obtained for the isothermal equation of state p = ke 
with m = 1, Q = 95/100 and A = —1/100. Dashed vertical lines denote locations of 
the horizons. 


and in particular 


(r - 1 - a^) 


\ 


2m A 


+ ^ + Koo) = 

’’oo 3 


= (r-1 - aDyi 


2 m Q2 
r r^ 6 


Using equation (ITTl) one can find, that for the polytropic equation of state 


P poo 


r-1 


1 

r -1 


r - 1 - a2 


(45) 


(46) 


This, in connection with equation flT^ . yields 


u = u. 


r_i 

oo / ^OO ^ ^ 


^ ^2 


Combining equations 

(r -1 - 4} 


Vo" r-i-a^; 

and fH7|) one gets 


\ 


2m Q2 A „ 

1 -h -r 2 + S = 

^ rL 3 ~ ^ 


= (r-i-oL) 


\ 


(47) 


3m Q'^ A 2 

' - 27 + ^ - I’-* + ■ 


(48) 


where 


B 





(49) 


Substituting al and with formulas taken from equation ([20]), one gets equation fjTOll 
with r* as the only one unknown. Solving this equation with boundary values of Too 
and ttoo gives positions of the critical points r* and, as a result, values of u*, al- The 
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knowledge of roo, floo, and al allows one to determine from equation fH7|) . 

Finally one can calculate the constant in equation (jH]), at roo- Combining this formula 
with equation (H7jl provides an implicit formula for u^{r). Solving it numerically and 
using equation (|8]) one can get the function (p//utYir). Similiary to the case of the 
isothermal fluid one may variate the constant B from equation obtaining solutions 
that are not transonic. Figure [5] shows a sample plot of the function {pf /utY{r). 

Figure [5] demonstrates, that it is possible to get a closed, subsonic solution. One 
can obtain a solution in a binocular shape consisting of two closed curves — one 
above the sonic point and one below it, laying partially under the event horizon. The 
transonic solution is also closed and similar to one obtained for isothermal fluid with 
k = \. It has been known, that in the case of the Reissner-Nordstrom spacetime 
one can obtain flow trajectories closed near the Cauchy and event horizons [H]. On 
the other hand, the possibility of getting trajectories closing for radii far from the 
event horizon in the Schwarzshild-anti-de Sitter spacetime has been revealed in [2]. It 
suggests, that in the presence of the both cosmological constant and charge it should be 
possible to get trajectories closed from both sides. The results described here conhrm 
that in special cases of the subrelativistic isothermal flow or polytropic flow. This is 
especially interesting in the context of rotating black holes. As it was mentioned in the 
Introduction, there exist qualitative similarities between Reissner-Nordstrom spacetime 
and the Kerr spacetime. For example one may apply swiss cheese model to a charged 
black hole in order to approximate real astrophysical situation [?]. The obtained results 
reveal interesting properties of the flows in such a situation. 

6. Summary 

In this paper I considered isothermal and polytropic flows in the Reissner-Nordstrom- 
(anti-)de Sitter spacetime. I described the way to hnd sonic points and to obtain the 
solutions {iB/utYir). One may think about these sonic points as critical points of a 
hamiltonian dynamical system connected with this flow. I investigated specihc cases 
of isothermal test fluids: so-called ultra-stiff fluid, ultrarelativistic fluid, radiation fluid, 
and subrelativistic fluid. In these cases the problem of solving the equations of motion of 
the fluid can be reduced to solving a polynomial of at most four degree. Exemplary plots 
of the flows were showed. I proved the correlation between positions of the horizons and 
the sonic points for the radiation fluid. Polytropic flows were also studied. It occured 
that in case of the subrelativistic isothermal flow or polytropic flow it is possible to obtain 
closed, binocular-like trajectories in (pf /utY{r) phase space. This property is consistent 
with previous results considering spacetimes with either non-zero cosmological constant 
[2] or charge [H], where trajectories were closing respectively in the neighbourhood of 
the cosmological horizon and the Cauchy horizon. 

This is especially interesting in the context of rotating black holes. As it was 
mentioned in the Introduction, there exist qualitative similiarities between Reissner- 
Nordstrom spacetime and the Kerr spacetime. The obtained results reveal interesting 
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Figure 5. Solutions for polytropic flow characterized by m = 1, Q = 9/10, 
A = —3.535 X 10“^, r = 4/3, = 10®, and = 0.2. The plot presents transonic 

solution (obtained for B = Bq) and two non-transonic. 


properties of the accretion in the Reissner-Nordstrom case. I believe that at least some 
of them can be present in Bondi-type accretion on rotating black holes. 

This paper fits into current research on closed flow trajectories (cf. [21 HI [TH |2T]). 
It appears, that whether solutions of this type exist, depends on both the spacetime 
and the equation of state of the fluid. In the case of the Reissner-Nordstrom-anti-de 
Sitter spacetime it is relatively easy to obtain such solutions — they appear to exist for 
matter models that can be interpreted as a gas of massive particles [1]. In other types of 
spacetimes it may be not so common, for example for a pure Schwarzschild solution one 
may have closed trajectories in more exotic case of polytropes with a high polytropic 
index [21]. The dependence of the existence of the closed solutions of considered 
spacetime and fluid types is a problem that needs further investigation. Of course 
it would be especially enlightening to compare results obtained here to corresponding 
results for Kerr spacetimes, which unfortunately do not exist so far. 

Acknowledgments 

Primarily I thank Dr. Patryk Mach for the assistance. Without his endless aid, long 
dicussions, and explanations, this paper would never come into existence. Also some 
numerical work was based on his code. I thank Jerzy Knopik for productive talks, which 












Bondi-type accretion in the Reissner-Nordstrdm-(anti-)de Sitter spacetime 


14 


sometimes were crucial for understanding. At the end I thank Prof. Edward Malec for 
important remarks. 


Appendix 


In this appendix I prove the relation between locations of the horizons and critical points 
that was discussed in Section 14.31 In the generic case of the Reissner-Nordstrom- 
de Sitter spacetime, when there exist three horizons, one of them is always located 
above both critical points, one is between the critical points, and one is below both 
critical points. Similiar property also holds for the Reissner-Nordstrom-anti-de Sitter 
spacetime: when there exist two horizons, one of them is between the critical points, 
and the second one is below both critical points. 

In order to prove it, let me consider two polynomials: 


fir) 

air) 


A 4 2 

= —r — r + 2 mr, 
3 


= —r H— mr. 


(A.l) 

(A. 2 ) 


These functions are formulas for derived from the equations ([2]) and ([37]). For hxed m 
and A, the function /(r) yields a value of for which there is a horizon in r. Similiary 
the function g{r) yields for which there is a critical point in r. Inverses of these 
functions were plotted in hgure O I consider behaviour of these polynomials for r > 0 
and > 0. The function /(r), as a quartic function, may have at most three extrema. 
The function g{r) is a quadratic function with a maximum at r = 3m/2. Sample plots 
of functions /(r) and g{r) are presented in hgure lATl As I demand the existence of three 
horizons in A > 0 case and two horizons in A < 0 case, I am interested in situations 
when for r > 0 and a hxed value of there are respectively three and two values of r 
such that /(r) = 

By diherentiating /(r) one can obtain a condition for extrema of /: 


tU ^ 3 

/ (^) = 


2 r + 2 m = 0 . 


(A.3) 


The case when an extremum is an inhection point if'fr) = 0) will be considered later. 
Substracting equations flA.ll and IA.2p results in 


fir)-air) 


r 

4 


/4A o 

( —r — 2 r + 2 m 


(A.4) 


It implies, that for r 7 ^ 0 functions /(r) and g{r) have /(r) = g{r) at extrema of fir). 
Both functions /(r) and g{r) vanish at r = 0. Moreover one has ffO) = 2m > ‘imj‘2 = 
5 ''( 0 ). Let me now consider the two cases. 

At hrst I will assume A > 0. I demand the existence of three horizons, so for r > 0 
the function /(r) has to have two extrema (it is clear that it cannot have more extrema, 
as /'(O) > 0). Beacause /'(O) > 5''(0) > 0 and A > 0 the hrst of these extrema will be a 
maximum (at ri) and the second one will be a minimum (at r< 2 ). Of course firf) = airi) 
and /(r 2 ) = air 2 )- It means that the function /(r) is increasing for r < ri. Then it is 
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decreasing until r reaches r 2 . After that /(r) goes to inhnity. Three horizons can exist 
only, if a hxed (such that for the horizons /(r) = Q"^) satishes r 2 < < ri. It 

implies, that one of the horizons is located on the left side of the parabole g{r) (below 
the critical points), one inside it (between the critical points), and one on the right side 
(above the critical points). 

A similiar argument works in the case of A < 0. This time I demand the existence 
of two horizons in the half-plane r > 0. Again one has f'{0) > g'{0) > 0. Beacause 
lim^-j-oo/(?") = —oo, the graph /(r) has to cross the graph g{r) in a point r such that 
0 < r < 3m/2. Then /(r) starts decreasing and cannot cross g{r) again. In the opposite 
situation it would have a minimum and then would start increasing, while g{r) would 
be still decreasing. The function /(r) could not change the sign of derivative and one 
would have limr_).oo f{r) = oo, what stays in a contradiction with the assumption that 
A < 0. It means, that there exists at most two horizons — one below both critical 
points and one between them. 

The function /(r) does not have an extremum in the point described by equation 
flA.3jl . if additionally /"(r) = 0 in the same point. Both of these conditions lead to 
A = 2/(9m^) and the occurence of a inflection point at r = 3m/2 (the maximum of 
the function g{r)). Obviously this is possible only in the Reissner-Nordstrom-de Sitter 
spacetime. In this case there is only one horizon (which is located below both critical 
points). 





Figure Al. Plots of /(r) (dashed) and g{r) (solid) for m = 1 and sample A values. 
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